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MINIMUM THRUST FOR CORRECTING KEPLERIAN ORBITS WITH 


APPLICATION TO INTERPLANETARY GUIDANCE 
By Frederick W. Boltz 
Ames Research Center 


SUMMARY 


An analytical investigation has been conducted to determine minimum 
thrust requirements for effecting small changes in several parameters of 
Keplerian orbits. Impulsive thrust and the approximations of linear perturba- 
tion theory have been assumed. The results are generally applicable to essen- 
tially two-body planar motion where the natural deviations of the vehicle 
flight path from a conic section are sufficiently small. 

The first-order thrust requirements obtained in the analysis are pre- 
sented in terms of the optimum thrust angle and the minimum derivative of 
thrust-produced velocity impulse with respect to the particular orbital param- 
eter to be corrected. Minimum thrust requirements are found for correction at 
a given point along an orbit of the semimajor axis, the eccentricity, and the 
true anomaly or orientation of the line of apsides (argument of pericenter) 
within the orbital plane. Minimum thrust requirements are also found for cor- 
recting the orbital radius at an arbitrary point and the elapsed flight time 
to arrive at that point. In addition, thrust requirements are obtained for 
correction at a given point along an orbit of the elapsed flight time to a 
second arbitrary point having a fixed value of Orbital radius. 

A general formulation of the thrust requirements for simultaneous correc- 
tion of any pair of the orbital parameters that define the in-plane motion is 
presented in terms of the minimum thrust requirements for correction of each 
of the parameters separately. Application of this result to the problem of 
implicit in-plane guidance during the midcourse and approach phases of inter- 
planetary flight is indicated. In particular, it is shown how the minimum- 
thrust results obtained in the present analysis may be applied to both fixed 
and variable time of arrival as well as fixed radius of pericenter types of 
guidance which make use of flight-path deviations from a reference trajectory. 


INTRODUCTION 


The problem of spacecraft guidance has been the subject of considerable 
research in recent years. This research has produced a number of methods for 
simplifying the onboard computation of the vehicle velocity correction. These 
methods may be generally classified as either implicit or explicit guidance, 
depending on whether or not a reference trajectory is used. Explicit guidance 
requires numerical integration of the equations of motion but has the advantage 
of greater flexibility in regard to off-design conditions. In implicit 



guidance, numerical integration is avoided by utilizing a reference trajec- 
tory. It is this type of guidance which is of primary interest in the present 
analysis. Corrections to the spacecraft velocity are calculated using pertur- 
bation expressions evaluated on the reference orbit in conjunction with pro- 
jected flight-path deviations from the reference orbit at the terminal point. 
The method is based on the assumptions that the velocity corrections are 
impulsive and that the variant motion about the reference trajectory is ade- 
quately described by linearized equations of motion. If two-body dynamics can 
be assumed, the required corrective velocity impulse can be obtained analyt- 
ical 1 y from parameter perturbation expressions of the type found in refer- 
ences 1 through 5. When only a single orbital parameter is to be corrected, 
the thrust angle can be selected to minimize the velocity impulse. The opti- 
mum use of impulsive thrust to correct the pericenter radius during approach 
to a planet has been considered in references 6 and 7. In general, it is pos- 
sible to correct three independent orbital parameters or three components of 
flight-path deviation from a target point with a single application of thrust. 
The problem of finding the required thrust impulse for this type of correction 
has been extensively treated in references 8 through 14. 

The present investigation is concerned with the derivation of minimum 
thrust requirements for single-parameter corrections and with the application 
of these results to several types of spacecraft guidance. A general formula- 
tion of minimum thrust requirements for correcting any orbital parameter has 
been obtained along with specific formulas for a number of parameters of inter- 
est. These results were obtained by applying the condition for minima given 
in the elementary calculus to general expressions for parameter perturbations 
due to impulsive velocity changes. This method is the same as that used in 
reference 7 for correcting pericenter radius and provides only first-order 
solutions. While the minimum- thrust requirements for correcting some param- 
eters (such as period of orbit and pericenter radius) are well known, this is 
not the case for other parameters considered in the analysis (such as orbital 
radius and flight time to an arbitrary point on the orbit) . 

With the admission of linearized two-body dynamics the general guidance 
correction can be separated into a correction of the plane of motion and an 
independent correction of motion within the orbital plane. The in-plane cor- 
rection can be generally treated as the correction of two independent orbital 
parameters. It is shown in the present analysis how the in-plane corrective 
velocity vector can be calculated from minimum thrust requirements for the 
various independent parameters. Several different types of midcourse and 
approach guidance for use during interplanetary flights considered include 
both fixed and variable time of arrival and fixed radius of pericenter (cf. 
refs. 11-14). 

Along with the assumptions of impulsive thrust and two-body dynamics, it 
is assumed in the analysis that the desired changes in the parameters to be 
controlled are sufficiently small to justify the use of linear perturbation 
theory. In most instances these assumptions provide sufficient accuracy in 
the computation of required velocity corrections and are almost a necessity 
when analytic solutions are desired. 
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NOTATION 


a semimajor axis of orbit, positive for elliptic orbits, negative for 

hyperbolic orbits 

b semiminor axis of orbit, positive for elliptic orbits, imaginary for 

hyperbolic orbits 

e eccentricity of orbit 

E eccentric anomaly, positive for ascent (0 <_ E <_ tt) , negative for 

descent (-tt ^ E 0) 

H hyperbolic anomaly, positive for ascent, negative for descent 

p semilatus rectum of conic section 

P general symbol for any orbital parameter of interest 

r radial distance from central body 

t time 

T period of orbit 

u horizontal component of velocity 

u c circular orbital velocity 

v vertical component of velocity 

V resultant velocity relative to inertial coordinates 

thrust-produced velocity impulse 

y flight-path angle relative to local horizontal, positive for ascent, 

negative for descent 

S( ) operator signifying small deviation 
A ( ) operator signifying increment or difference 

6 true anomaly, positive for ascent (0 0 tt) , negative for descent 

(-7T < 0 < 0) 

p product of universal constant of gravitation and mass of celestial 

body 

x thrust angle relative to local horizontal, positive for thrust vector 

increasing vertical component of velocity, negative for thrust 
vector decreasing vertical component of velocity 
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T 


a 


T 


T 

e 


T 


T 


0 


T 


e 0 


T 


r 


T 


T 


I 

r 


thrust angle to maximize change in semimajor axis or period of orbit 

thrust angle to minimize change in semimajor axis or period of orbit 

thrust angle to maximize change in eccentricity 

thrust angle to minimize change in eccentricity 

thrust angle to maximize change in true anomaly 

thrust angle to minimize change in true anomaly 

thrust angle to maximize change in orbital radius at fixed central 
angle from thrusting point 

thrust angle to minimize change in orbital radius at fixed central 
angle from thrusting point; also, thrust angle to correct flight 
time to fixed point on orbit 

thrust angle to maximize change in orbital radius having fixed value 
of true anomaly 


thrust angle to maximize change in flight time to traverse fixed 
central angle 

x t thrust angle to minimize change in flight time to traverse fixed 

° central angle 

t » thrust angle to maximize change in elapsed flight time to arrive at 
fixed value of true anomaly 

(j) angle from pericenter to any point on orbit measured at center of 

orbit 


( ) quantity normalized with local circular orbital velocity, /jj/r 


Subscripts 

a value at apocenter 

p value at pericenter 

1 value at orbital position where thrust impulse is to be applied 

2 value at orbital position where primary change in orbit is to be 

effected 


4 



ANALYSIS 


Assumptions and Restrictions 

It is assumed in the analysis that vehicle thrust produces an impulsive 
velocity increment A*V X which is added vectorially to the velocity vector of 
the vehicle. It is also assumed that small perturbations in the vehicle 
motion due to impulsive thrust may be calculated with sufficient accuracy by 
means of linearized two-body, point-mass equations of motion. The restric- 
tions on orbital conditions in the case of a noncentral force field require 
simply that the perturbations of the vehicle flight path from a conic section 
be relatively small. Within the limitation of these restrictions, the linear- 
ized two-body orbital relations may also be used as a basis for determining 
minimum thrust requirements for effecting small changes in certain orbital 
elements or trajectory variables. 

One effect of the restriction in size of impulsive thrust is that pertur- 
bations of the vehicle motion within and normal to the orbital plane are inde- 
pendent. Since thrust directed normal to the orbital plane does not affect 
motion within the plane of orbit (to first order) , the present analysis is 
restricted to in-plane thrusting only and the resulting effects on the in-plane 
motion. 

The convention of notation used in the present analysis is illustrated in 
figure 1, which indicates the geometric parameters of elliptic and hyperbolic 
orbits pertinent to the analysis. Parabolic orbits may be considered a spe- 
cial case of elliptic orbits in which the semimajor axis, a, is infinitely 
large. For hyperbolic orbits a is negative but has a positive counterpart 
in the semi transverse axis, -a. 

In the derivation of analytical expressions that follows, intermediate 
steps are often omitted or indicated only briefly. Since in many instances 
the auxiliary relations used in the derivation are not given, a short summary 
of two-body point-mass relations pertinent to the present analysis is presented 
in appendix A. (A more detailed summary is to be found in ref. 2.) 


General Formulas for Optimum Correction of Orbital Parameters 

Impulsive thrust applied to an orbiting vehicle produces an instantaneous 
change in the local velocity vector and related changes in all orbit param- 
eters. If P is any parameter of interest, it can be expressed as a function 
of the local velocity vector and the orbital radius at the point of thrust 
application. Thus, with the velocity vector described by its magnitude and 
direction (V ls Yi) or by its horizontal and vertical components (U} , v^). 


p = f(Vi, Yi, Ti) = g(uj, Vi, ri) 


( 1 ) 
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The perturbation of P due to a small change in the local velocity vector is 
given by 


dp = wr dVl + wr dY1 = l£r dui + dVi 


( 2 ) 


since the orbital radius at the thrusting point is considered fixed. 



Sketch (a} 


The changes in vehicle velocity 
due to the application of impulsive 
thrust are illustrated in sketch (a) . 
The velocity increment associated with 
the thrust impulse is described by its 
magnitude, AV T , and direction, t, 
referred to the local horizontal. It 
is readily seen that, in the limit as 
AV t goes to 0, changes in the various 
components of the velocity vector are 
given by 


d\l i = dV T cos(x - yi) 

dv T 

dYi = sin(T - Yi) 

duj = dV T cos t 
dvi = dV T sin x 


J 


(3) 


where x - yi is the thrust angle relative to the flight path or velocity 
vector. When these differential expressions are substituted into equation (2), 
it is found that 


or 


dP 


3f 

_3Vi 


cos (x 


Yl) 


+ 


1_ 9f_ 
Vi 9 yi 


sin(x 



dV 

x 


(4a) 


dP 


9g_ 

3ui 


cos x 


+ 


9g_ 

9v 2 


sin x dV 


c 


(4b) 


These equations provide alternate expressions for the derivative dV T /dP in 
terms of thrust angle relative to the flight path or relative to the horizontal. 
This derivative may be used to obtain a first-order estimate of the thrust mag- 
nitude, AV t , required to change the parameter P an amount AP according to 
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AV 


T 



AP 


The thrust magnitude will be a minimum when the thrust angle, x, is selected 
so as to minimize the derivative. This optimum thrust angle is found by set- 
ting the derivative of equation (4a) or (4b) with respect to x equal to 0 
and solving for x. If x* represents the optimum thrust angle, it is found 
that 


or 


tan(x* - Yl ) = 


tan x* = 


9f /9y i 

(5a) 

Vi af/aVi 

9g/9Vi 

9g/9u! 

(5b) 


from which two values of x* - and of x* (differing by 180°) are obtained. 

The corresponding minimum value of dV /dP is obtained by substituting 

L * 

x* for x in equation (4a) or (4b) with the result that 


or 



(6a) 


(6b) 


where the plus and minus signs correspond to thrusting in opposite directions. 

The thrust angle for which the parameter P is invariant is found by 
setting equation (4a) or (4b) equal to 0 and solving for x. If x*^ 

represents this thrust angle, it is found that 


tan x* 

o 


3g/3Uj 

9g/9Vi " _COt T * 


(7) 


so that, to first order, the thrust angle for no change in a given parameter 
is different by 90° from that for maximum change. It should be noted that 
for thrust at this angle the parameter does in fact change, but by an amount 
proportional to second- or higher-order derivatives. 

It follows from the above results that the general linearized expression 
for the parameter perturbation (eq. (4)) can also be given, in terms of the 
optimum thrust angle and the maximum derivative of the parameter with respect 
to the thrust-produced velocity impulse, by 


7 



dp 


dV T cos(t - t*) 


( 8 ) 


dP 


dV„ 


since, to first order, only the component of thrust in the optimum- thrust 
direction contributes to the perturbation. 

The first part of the following analysis is concerned with developing 
linear perturbation expressions for various parameters of interest. In the 
second part minimum thrust requirements are determined for the control of 
these parameters. 


Perturbation Formulas for Various Orbital Parameters 

The perturbative effects of impulsive thrust on various elements or 
parameters of an orbit are obtained from the general linearized perturbation 
expressions of the previous section and the appropriate two-body relations for 
these parameters. The orbit elements considered in the present analysis are 
the semimajor axis (or period of orbit), the eccentricity, and the true anom- 
aly at the point of thrust application. Perturbations of the orbital radius 
at any other point on the orbit and of the elapsed flight time to any other 
point are also considered. Perturbations of the semimajor axis and the eccen- 
tricity reflect changes in the orbit size and shape. The perturbation of true 
anomaly at a given point is equal to the rotation of the line of apsides of 
the orbit within the orbital plane. 

In the following presentation of perturbation expressions, alternate 
forms of the results are given in the case of semimajor axis (or period of 
orbit), eccentricity, and true anomaly. These forms differ in the specifica- 
tion of thrust angle (either measured with respect to the flight path or the 
horizontal) and in the particular component of velocity with which the thrust- 
produced velocity differential is normalized. For perturbations of orbital 
radius and elapsed flight time, only one form is presented here; the other in 
each case is given in appendixes B and C, respectively. 

Perturbation of semimajor axis or period of orbit .- The perturbative 
effects of impulsive thrust at orbital position (r la 0 ^ ) on the semimajor axis, 
a, or the period of orbit, T, are obtained by applying equation (4) to the 
two-body relations for these elements given in appendix A. It is found that 

Aa 9 AT 9 ^ dV T 

— = 3 ■ Y~ = 7 — ~^2 [(1 + 2e cos 0! + e 2 ) cos (t - Yl)] ( 9a ) 

dT 2(1 + e cos 0i) dV T 

=— = 5 [(1 + e cos 0i)cos t + e sin 0* sin t] 

1 1 - e 2 U 1 

(9b) 

The value of a is positive for elliptic orbits and negative for hyperbolic 
orbits. Both a and da are undefined for parabolic orbits. 


or 

da 

a 
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Perturbation of eccentricity .- The perturbative effect of impulsive 
thrust at orbital position (ri, 0 i ) on the eccentricity of orbit, e, is 
obtained by applying equation (4) to the two-body relations for this element 
given in appendix A. It is found that 

f (1 - e 2 ) sin 0! I dV T 

de = [_2Ce + cos 0i)cos(t - Yl ) + : + e cos sin(r - Yi)J ^ 

(10a) 


or 

dV T 

de = {[(2 + e cos 0i)cos Gj + e]cos t + (1 + e cos 0i)sin 0! sin x} 


(10b) 

In the case of nearly circular orbits (e -* 0) equations (10a) and (10b) 
reduce to 


de = (2 cos 0i cos x + sin 0i sin x)dV 

x 


(10c) 


This expression cannot be used directly for perfectly circular orbits, since 
the true anomaly is undefined. However, since the change of eccentricity is 
the value of eccentricity after thrusting, it is given by 

de = A - (2 - V 2 )!! 2 (lOd) 


Here V and u are the normalized values of velocity and horizontal component 
of velocity after thrusting which, with Vx = u\ = 1 and yi = 0 (see 
sketch (a)), are given by 

V = A + (dv x ) 2 + 2dV T cos T 
u = 1 + dV T cos x 


Thus, for a circular orbit it is found that, to first order. 


de 


- /T7 


3 cos 2 x dV 


(10e) 


It is interesting to note that this result can also be obtained from 
equation (10c), if the value of true anomaly attained at the thrusting^ point 
after the application of thrust is substituted for 0j. With V and u as 


9 



given above, this value of true anomaly is readily obtained from the two-body 
relation 


cos 0 = 

It is found that, to first order, 

cos 0 


u* 


/l - (2 - V 2 )u^ 


2 COS T 


A + 3 


COS 2 X 


with dependence only on thrust angle. 


Perturbation of true anomaly .- The perturbative effect of impulsive 
thrust at orbital position (ri, 0i) on the true anomaly, 9^, is obtained by 
applying equation (4) to the two-body relations for this parameter given in 
appendix A. It is found that 


d0 l " - ? 


2 sin 0i cos(x - yi) - 


2e + (1 + e 2 ) cos 0i 
1 + e cos 0i 


sin(x - yi) 


dV„ 


V i 


(11a) 


or 


dV 


1 T 

d0i = - — [(2 + e cos 0i)sin 0i cos x - (1 + e cos 0i)cos 0i sin x] 

e U i 


(Hb) 


This change of true anomaly at the thrusting point is equal to the change of 
true anomaly at any other point, since it is due to a slight rotation of the 
line of apsides of the orbit within the orbital plane. 

Perturbation of orbital radius .- The perturbative effect of impulsive 
thrust at orbital position (ri, 0 i ) on the orbital radius at position 
( r 2 , e 2) is obtained from the conic relation 

1 + e cos 0 a 
1 + e cos 0 2 / 

The logarithmic derivative of this expression, with ri constant, gives 

e sin 0! d0 i e sin 0 2 d0 2 

1 + e cos 0i 1 + e cos 0 2 


dr 5 


cos 0! 


cos 0 2 


1 + e cos 0i 1 + e cos 0; 


de 



( 12 ) 
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For ri fixed, the differentials de and dQ\ are given by equations (10) and 
(11), respectively, in terms of e, 0i, t, and dV T /V 1 . The differential d0 2 
can be set equal to d0i or 0 depending on the particular problem under con- 
sideration. If the central angle 0 2 - 0i is considered fixed, then d0 2 
is equal to d0i as given by equation (11). In this case equation (12) gives 
the perturbation of the orbital radius r 2 along an inertial ly fixed radial 
line in the orbital plane. If the value of 0 2 is considered fixed, then 
de 2 is equal to 0, and equation (12) gives the perturbation of the orbital 
radius r 2 along a line which, in general, undergoes rotation with thrust so 
as to maintain a fixed value of true anomaly, 0 2 . Since the coefficient of 
d0 2 in equation (12) vanishes for 0 2 equal to 180° or 0, perturbations of 
apocenter and pericenter radii can be determined without specification of d0 2 . 

For fixed central angle, 0 2 - 0^, equation (12) can be expressed as 


dr 2 = 1 
r 2 1 + e cos 0 2 


{[2(1 - cos A0 ) - e sin A0 sin e^cos t 


+ (1 + e cos 0i)sin A0 sin t} 


^2 

Ul 


(13a) 


and, for fixed true anomaly, 0 2 , as 
dr„ t 


*2 


1 + e cos Q' 


{[2(1 - cos 0i cos 0 2 ) + e sin 2 0j cos 02 ]cos t 


+ (1 + e cos 0i)sin 0! cos 0 2 } 


£ 4 

u i 


(13b) 


Alternate forms of these perturbation expressions in terms of thrust angle 
relative to the flight path, x - y ly are given in appendix B. 

Perturbation of elapsed flight time .- The perturbative effect of impul- 
sive thrust at orbital position (r 1? on the elapsed time of flight to 

position (r 2 , 0 2 ) is obtained from Kepler’s equation expressed as 


t 2 ~ t i =/— [ E 2 - E 1 - e ( sin e 2 " sin E x )] 


where E is the eccentric anomaly. The logarithmic derivative of this 
equation gives 


d(t 2 " t l) 
1 2 - t l 


3 da t /a 3 /y r f1 
2 a t 2 - t x LU 

- (1 “ e cos E|)dE^ 


e cos E 2 )dE 2 
(sin E 2 - sin E^de] 


(14) 
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The differentials da and de are given, for r x fixed, by equations (9) and 
(10), respectively, in terms of e, 0 X , t, and dV T /V x . The remaining differen- 
tials, dE 2 and dE x , are obtained from formulas given in appendix A as 



1 


A - e2 

+ e cos 


e. 

j 



de 

1 - e2 


+ de 


3 > 


j = 1,2 


where de and d0 x are given by equations (10) and (11), respectively. The 
single unspecified differential, d0 2 , can be set equal to d0 x or to 0 
depending on whether 0 2 - 0 X or 0 2 is considered fixed. In the former case 
the terminal point remains inertially fixed, and in the latter case it does 
not. Both constraints are useful in guidance computations as is subsequently 
shown . 


For fixed central angle, 0 2 - 0 X , equation (14) can be expressed as 


d(t 2 - t x ) 


A» 3 /m 

to - tl 


1 + e cos 0i 


X x cos x + X 2 sin t 


dV. 


Ul 


(15) 


where 


*1 

x 2 


A - B(2 + e cos 0 x )sin 0 X 
1 + e cos 0 X 


+ K ( 1 + e cos 0 X ) 


B cos 0 X + Ke sin 0 X 


with 

(2 + e cos 0 2 )sin 0 2 (2 + e cos 0 x )sin 0 X 

A 

(1 + e cos 0 2 ) 2 (1 + e cos 0 X ) 2 

B = (1 + e cos 0 2 ) -2 - (1 + e cos 0 X ) -2 



For fixed true anomaly, 0 2 , B is replaced in the above expressions by 


B' = -(1 + e cos 0 1 ) 2 


These results are valid for both elliptic and hyperbolic orbits, since 
Kepler's equation is applicable to any type of orbit. It is only the form of 
this equation that need be changed for use in the different cases. This is 
simply done by using the identity E = -iH or H = iE (cf . appendix A) . 
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However, the above results, which do not contain any imaginary terms for cer- 
tain values of eccentricity, are independent of the usual restriction on type 
of orbit specified. 

In the case of parabolic orbits, the factor K in equation (15) is 
indeterminate, and a limiting value must be obtained. If this limiting value 
is represented by L, it is found that 


L = lim K 
e+l 


0 2 

0i 

i 

/ 5 0 2 

5 01 YI 

tan — - 

tan — - 

5 

(tan 5 — - 

tan 2~ )_ 


In the case of circular orbits, equation (15) reduces to the limiting 

value 


d(t 2 - tj) 
t 2 - t\ 


= tt [3 A0 - 4 sin A0)cos 
AU 


t + 2(1 - cos A6)sin x] 


dV 

T 


Alternate forms of these perturbation expressions in terms of thrust 
angle relative to the flight path, x - yi, are given in appendix C. 


Minimum Thrust Requirements for Correction of Various 

Orbital Parameters 

The minimum thrust requirement as used throughout this report refers to 
the optimum thrust angle for correcting a particular orbital parameter and the 
minimum derivative of thrust-produced velocity impulse with respect to that 
parameter. Such minimum thrust requirements for correcting various parameters 
of an orbit are obtained by applying equations (5) and (6) to the appropriate 
two-body relations for these parameters. The partial derivatives appearing in 
equations (5) and (6) may be obtained from a comparison of equation (4) with 
the perturbation expressions for the various parameters given in the previous 
section. 

In general, there is a particular point along an orbit where the minimum 

derivative, dV T /dP, attains an absolute minimum value. This optimum orbital 
* 

location for parameter correction is found by setting the derivative of 
dV /dP with respect to 0i equal to 0 and solving for 0j. Such optimum 

values of 0^ are also presented for several of the orbital parameters 
considered in the analysis. 

As in the previous section, alternate forms of the results are given in 
the case of semimajor axis (or period of orbit), eccentricity, and true anom- 
aly. In the case of orbital radius and elapsed flight time, only one form is 
presented here; the other in each case is given in appendixes B and C, 
respectively . 

Correc tion of semimajor axis or period of orbit .- If the thrust angle for 
maximum change in the semimajor axis (or period of orbit), for a given small 
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(■nil 


amount of thrust, is represented by x a , it is found from equation (5) that 


tan(T a - yi) = 0 


(16a) 


or 


tan t & = tan yi 


e sin 0j 
1 + e cos 0i 


(16b) 


It is readily seen that this optimum thrust direction is valid regardless of 
the amount of thrust applied, since in all cases a maximum change in Vj is 
desired for a given value of dV T . No change in semimajor axis or period 
occurs if the speed, Vi, is unchanged by thrusting; this requires thrust 
normal to the flight path or velocity vector for small amounts of thrust. 


The minimum dimensionless 1 derivative of thrust-produced velocity impulse 
with respect to semimajor axis or period of orbit is found by substituting 
x a for t in equation (9) to obtain 


a 

Vi 



3 T dVx a 
2 Vi dT 


+ 


1 - e 2 

2(1 + 2e cos 0i + e 2 ) 


(17a) 


or 


a 




3 T dV -a 
2 V p dT 


1 - e 


2v/f + 2e cos 0i + e 2 


(17b) 


From equation (17b) it is found that the absolute minimum values of 
dV T /da and dV Ta /dT occur at a value of 0i equal to 0. 

Correction of eccentricity .- If the thrust angle for maximum change in 
the eccentricity of orbit, for a given small amount of thrust, is represented 
by r e , it is found from equation (5) that 


(1 - e 2 )sin 0i 

tan( T e 7i) 2(1 + e cos 0i)(e + cos 0^) 


(18a) 


or 


(1 + e cos 0i)sin 0i 
e (2 + e cos 0i)cos 0j + e 


(18b) 


^he differential dV T can be normalized with any convenient speed param- 
eter. In this report a variety of such parameters are used for this purpose 
including Vi, V 2 , Ui, V c , etc. 
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The thrust angle for which the eccentricity is invariant, to first order, is 
different from x e by 90°. 

The minimum dimensionless derivative of thrust-produced velocity impulse 
with respect to eccentricity is found by substituting x e for x in equa- 
tion (10) to obtain 

, dV T r r(l - e 2 )sin 0 i"1 2 T 1//2 

VT^ = 1 { 4(e + cos 01)2 + [ l + e cos eT J } (19a) 

or 

1 dV T P , o O - 1 / 2 

— — = ±{[(2 + e cos 0j)cos 0i + e] 2 + (1 + e cos 0i) 2 sin 2 0i> 

(19b) 

For circular orbits, it can be shown from equation (lOe) that x e is 
equal to 0 or 180°; consequently, 


]_^Ie fie 

Vi de de 


±0.5 


(19c) 


From equation (19b) it is found that the absolute minimum value of 
dV T /de occurs at ©i equal to 0 or 180° with the absolute maximum value 

occurring at a value of ©i given by 


e 2 cos 3 0i + 3e cos 2 0i + 3 cos 0i + e(2 - e 2 ) = 0 


The one real root of this cubic equation for positive values of eccentricity 
is 

cos 0i = [(1 - e 2 ) 2 / 3 - 1] 

Correction of true anomaly .- If the thrust angle for maximum change in 
the true anomaly (i.e., maximum rotation of the line of apsides within the 
orbital plane), for a given small amount of thrust, is represented by x ft , it 
is found from equation (5) that 


tan (xq - yi) 


2e + (1 + e 2 )cos ©i 
2(1 + e cos ,0i)sin 0 


(20a) 


or 
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(1 + e cos 8i)cos 0i 
(2 + e cos 6i)sin Bj 


(20b) 


tan t 


The thrust angle for which the true anomaly is invariant, to first order, is 

different from t by 90°. 

0 


The minimum dimensionless derivative of thrust-produced velocity impulse 
with respect to true anomaly is found by substituting for t in 

equation (11) to obtain 


1_ 

Vi 



±e ^4 sin 2 0j_ + 


’2e + (1 + e 2 )cos 0i 
1 + e cos 0i 


2 -l/2 


(21a) 


or 



Ui d0 1 


±e[(3 + 2e cos 0i)sin 2 0^ + (1 + e cos © 1 ) 2 ] 


1/2 


(21b) 


From equation (21b) it is found that the absolute minimum value of 
dV T /d0i occurs at a value of 0j given by 

L 0 


e 2 cos 3 0i + 3e cos 2 0i + (3 + e 2 )cos 0i + 2e = 0 


The one real root of this cubic equation for positive values of eccentricity 
is 


COS 01= 




+ 




1 


Correction of orbital radius .- The thrust angles resulting in minimum 
required thrust for making a correction at a given point (ri, 0i) along an 
orbit of the orbital radius at a second point (r 2 , 0 2 ) are denoted by x r and 
x r , depending on whether 02 - 0i or 02 is considered fixed. These opti- 
mum thrust angles and corresponding minimum dimensionless derivatives of thrust- 
produced velocity impulse with respect to orbital radius are obtained from 
equations (5) and (6) with reference to equations (4) and (13) . For fixed 
central angle, ©2 - ©i, it is found that 

(1 + e cos ©i)sin A0 

tan r " 2(1 - cos A0) - e sin A0 sin 0i ^2) 


and 
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r 2 dV T 
T r 

ui dr 2 


1 + e cos 0 ; 


= ± 


A 


2(1 - cos A0) - e sin A0 sin 0 1 ] ^ + (1 + e cos 6 i ) 2 sin 2 A0 


(23a) 


or 


r 2 dV n 


r. _ 


1 + e cos 0 i 


U 2 dr2 




A 


2(1 - cos A0) - e sin A0 sin 0]J 2 + (1 + e cos 0i ) 2 sin 2 A0 


(23b) 


In the case of fixed true anomaly, 02 , it is found that 

(1 + e cos 0 i)sin 0 i cos 02 


tan x r ' = 


2(1 - cos 0 i cos © 2 ) + e sin 2 0 i cos ©2 


(24) 


and 


r 2 dV Tr , 

Ui dr 2 


1 + e cos 62 


= + 


^ 2(1 - cos 0 i cos 0 2 ) + e sin 2 0 i cos 0 2 ] 2 + (1 + e cos 0 i ) 2 sin 2 0 ] cos 2 0 2 

(25a) 


or 


r z dV T , 
L r 

u 2 dr 2 


1 + e cos 0 i 


= + - 


/[ 2(1 - cos 0 i cos 0 2 ) + e sin 2 0 i cos 0 2 ] 2 + (1 + e cos 0 1 ) 2 sin 2 0 X cos 2 0 2 

(25b) 


Alternate formulations of these minimum thrust requirements for orbital 
radius correction are given in appendix B. Formulas for the special cases of 
apocenter-radius and pericenter-radius correction are listed in appendix D. 
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Correction of elapsed flight time .- The thrust angles resulting in mini- 
mum required thrust for making a correction at a given point (r^ along 

an orbit of the elapsed flight time to a second point (r 2 , 0 2 ) are denoted by 
and T-t-f, depending on whether 0 2 - or 0 2 is considered fixed. These 

optimum thrust angles and corresponding minimum dimensionless derivatives of 
thrust-produced velocity impulse with respect to elapsed flight time are 
obtained from equations (5) and (6) with reference to equations (4) and (15) . 
For fixed central angle, 0 2 - 0i, it is found that 


and 


-L 1 

tan Tt = XT = A - B(2 + e cos e^sin 0! 

z + K(1 + e cos 0] 

1 + e cos 0! v J 


t 2 - tx dV Tt 
u i d(t 2 - tT) 


+ 


t 2 - ti e/(l + e cos ©j) 

j/frVu + X 2 


(26) 


(27) 


where Xj, X 2 , A, B, and K are functions of e, 0j, and 0 2 previously 
defined with equation (15). In the case of fixed true anomaly, 0 2 , the opti- 
mum thrust requirements are the same as those above except that B is 
replaced by B' . For parabolic orbits K is indeterminate and is replaced 
by its limiting value, L, given with equation (15). For circular orbits 
equations (26) and (27) in the case of fixed central angle reduce to the 
limiting forms 


tan t = 2 - (1 ~ co - 

1 t 3 A0 - 4 sin A0 


and 


t 2 “ "t i dV T j j 

— r? t? £— r- = ±A0//(3 A0 - 4 sin A0) 2 + 4(1 - cos A0) 2 

V l d(t 2 - ti) 

with dependence only on the central angle, A0. The thrust requirements for 
correction of elapsed flight time with fixed true anomaly at the terminal 
point are ambiguous for exactly circular orbits, since the true anomaly is 
undefined. 

In general, correction of elapsed flight time with the above minimum 
thrust specifications will result in a change in the orbital radius at the 
terminal point. However, there is one direction in which thrust will not 
change the terminal radius (to first order) . This thrust direction is normal 
to that required for maximizing the change in orbital radius. Thus, the 
thrust angle for correction of elapsed flight time with invariant terminal 
radius is either x r or x r ^, depending on whether 0 2 - 0j or 0 2 is 
considered fixed. T^ese thrust angles can be expressed as 
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tan T rQ = -cot t r (28a) 

tan x ! = -cot x 1 (28b) 

1 o 


where x^ and x r ’ are given by equations (22) and (24), respectively. 

When the thrust angles T r and x r * are substituted for x in 
equation (15), it is found that? for fixed central angle, 0 2 - 0^, 


t 2 - ti 


dV„ 


t 2 " 


e /sin 2 A0 + 


r2Cl - cos A0) - e sin A0 sin Qi"] 2 
|_ 1 + e cos 0 x J 


ui d(t 2 - t x ) /p 3 /P A sin A0 + 2B(cos 0 2 - cos ©x) + K[2e(sin 0 2 - sin 0x) + (1 + e 2 )sin A0] 

(29a) 

and, for fixed true anomaly, 0 2 , 


t 2 - tx 


dV T , 
r o 


t 2 - tx 


e /sin 2 0 x cos 2 0 2 + 


r2 (1 - cos 0 x cos 0 2 ) + e sin 2 0x cos 0 2 1 2 


1 + e cos 0x 


ux d(t 2 - t x ) /p 3 /P A sin 0 x cos 0 2 - 2B' (cos 0 2 - cos © x ) + K[2e + (1 + e 2 )cos 0 2 ]sin 0x 


(29b) 

where, again. A, B, B 1 , and K are functions of e, 0j, and 0 2 previously 
defined, with K replaced by its limiting value, L, in the case of parabolic 
orbits. For circular orbits equation (29a) reduces to the limiting form 



Vi d(t 2 - t^ 


A0>/sin 


2 


A0 + 4(1 - cos AG) 2 


3 A0 sin A0 - 8(1 - cos A0) 


with dependence only on the central angle, A0. The derivative expressed by 
equation (29b) is 0 for perfectly circular orbits but is meaningless in this 
case due to the ambiguity of the true anomaly. 

Alternate formulations of these thrust requirements for correction of 
elapsed flight time are given in appendix C. Formulas for the special cases 
of correction of flight time to apocenter and pericenter are listed in 
appendix E. 


Thrust Requirements for Simultaneous Correction of 
Two Orbital Parameters 

In general, it is possible to correct two independent orbital elements 
or parameters with a single application of impulsive thrust, since two inde- 
pendent scalar quantities are required to define the corrective thrust vector 
in the orbital plane. The specification of the corrective thrust vector in 
terms of the thrust angle, x, and the thrust magnitude, AV T , is obtained from 
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the simultaneous solution of the perturbation expressions for the two param- 
eters. Two different formulations for t and AV T based on the use of 
equations (4) and (8) are developed below. 

General formulation of thrust requirements .- Consider two independent 
orbital parameters. Pi and P 2 , requiring corrections 6P1 and <SP 2 . From 
equation (4b) it follows that 


&? 1 6P 2 

A v t = 

(Sgl/9u s )cos t + (9g 1 /3v s )sin x (9g 2 /3u s )cos x + (9g 2 /9v s )sin x 


where u g and v s are the horizontal and vertical (radial) components of the 
vehicle velocity at the time of thrusting. When these two equations are 
solved for x and AV X , it is found that 


tan x 


3g2 

9u~ 


6Pi 


3gl 

9v_ 


<5Pc 


3 gl 

3u_ 


<5P: 


3 g2 

9v„ 


6Pi 


(30) 


and 


AV, 



9 gl 9 g2 9 g2 9 gl 

9117 9v7 ' 9^7 9vJ 


(31) 


It is to be noted that two values of thrust angle differing by 180° are 
obtained from equation (30) . The value of x corresponding to a positive 
value of AV t is that for which the sign of (9gj/9u s )cos x + (3gj/9v s )sin x 
is equal to the sign of 6P • . The condition for the two orbital parameters 
being independent is that the denominator in equation (31) is not 0. In gen- 
eral, two parameters which are usually independent can become dependent at one 
or more points on an orbit. These points are specified by the conditions 
under which the denominator of equation (31) is 0. 


An alternate formulation of the thrust requirement for simultaneous 
correction of two independent orbital parameters may be obtained in terms of 
optimum thrust angles and minimum derivatives of thrust-produced velocity 
impulse with respect to each of the parameters. From equation (8) it follows 
that 


AV 


x 



6 P ! 



<SP 2 


dPj cos (x - x.^) dP 2 cos (x - t* 2 ) 
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where x* and x* are the optimum thrust angles and dV /dPi and 

12 * i 

dV /dP£ are the minimum derivatives. When these two equations are solved 
*2 

for x and AV^, it is found that 


AV_ cos x* - AV. cos x* 

T * x *2 t * 2 *1 

tan x = 

AV t sin x. - AV. sin x. 

t * 2 1 T *l 2 


(32) 


and 


AV„ 


AV = 
x 


,) 2 * (%) 2 


2 AV. AV cos (x ^ - x* ) 

T * 1 t * 2 2 *1 


sin(x - x ) 


(33) 


where 


dV 


AV. 


' * . 

j 


dP 


-i- 6Pj, j = 1, 2 


This thrust requirement can be expressed more concisely as 


AV X - AV X cos (x *2 - T * x ) 


tan (x - x*,) = 


AV sin (x - x ) 

t j. a a Tr 1 


(34) 


and 


T *i 

AV x = i (35) 

cos (x - X* x ) 

As before, two values of x differing by 180° are obtained from either equa- 
tion (32) or equation (34) . The value of x corresponding to a positive 
value of AV t is that for which |x - x #1 | or |x - x* | is less than 90°. 

The condition for the two orbital parameters being independent is that the 
optimum thrust directions for correction of the two parameters separately do 
not coincide. It is also to be noted that when the required correction in one 
of the parameters is 0, the thrust specification is in the direction which, to 
first order, does not change this parameter. Thus, the thrust requirement in 
this case is different from the minimum thrust requirement for correction of 
the other parameter singly. 
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Geometric interpretation of results . - A geometric interpretation of the 
results obtained above for computing the single velocity impulse required to 
correct two different orbital parameters using their minimum thrust require- 
ments is obtained from considera- 


tion of the vector diagram 



presented in sketch (b) . As is 
indicated in this sketch, the 
required velocity impulse may be 
resolved into either of two par- 
ticular pairs of orthogonal com- 
ponents. One of the two components 
in each pair is in the thrust 
direction required for maximum 
change of one of the two param- 
eters being corrected. The other 
component in each pair, being 
orthogonal to the first, has no 
effect, to first order, on the 
pertinent parameter. Thus, it is 
clear that, when the two thrust 
directions for maximum thrust 
effectiveness in the case of the 
two parameters are specified 
along with the velocity impulses 
AV and AV required in these 
*1 *2 

respective directions to correct 


each parameter separately, the 


velocity vector for simultaneous correction of the two parameters is uniquely 
defined. 


From sketch (b) it is apparent that the thrust vector for simultaneous 
correction of two parameters can be easily obtained graphically from the mini- 
mum thrust vectors for correction of the two parameters separately. All that 
is required is the erection of perpendiculars to the ends of these minimum 
thrust vectors. The intersection of these perpendiculars is the terminus of 
the desired single thrust vector correcting both parameters simultaneously. 

Also indicated in the vector diagram of sketch (b) is the geometric 
basis for the relations expressed in equations (34) and (35). It can easily 
be shown that, as is indicated, the angle t - x* between the vectors of 

magnitude AV X and AV T is the same as the angle between the line joining 

*1 

the ends of the vectors of magnitude AV and AV_ and the perpendicular 

T *l ~ *2 

to the vector of magnitude AV . With this correspondence in angles the 

X *2 

relation given by equation (34) is readily obtained. 
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RESULTS AND DISCUSSION 


Graphical Presentation of Results 

Minimum thrust requirements for correcting various orbital parameters 
are presented in figures 2 through 8. These results are given in dimension- 
less form and can be used for analyzing thrust -perturbed orbital motion 
around any celestial body. The formulas comprising the minimum thrust require- 
ments (optimum thrust angle and minimum derivative of thrust -produced velocity 
impulse with respect to the parameter of interest) are listed in table I for 
convenience. These thrust specifications are generally expressed as functions 
of the eccentricity of orbit and the true anomaly at the thrusting point. 

Values of eccentricity from 0 to ® have been considered as well as the full 

range of positive values of true anomaly. 

Optimum correction of semimajor axis or period .- The minimum thrust 
requirements for correcting the semimajor axis of various orbits are presented 
in figure 2. They are proportional to those for correcting the period of 
orbit. The optimum thrust angle for increasing both the semimajor axis and 

the period is equal to the flight-path angle. Thrust in the opposite direc- 

tion maximizes the reduction of these parameters. The variation of the flight- 
path angle with true anomaly is presented in figure 2(a) for various values 
of eccentricity. The extreme values of flight-path angle for elliptic orbits 
occur at the two midpoints of orbit (where the orbital radius is equal to the 
semimajor axis and the local velocity is equal to the local circular speed). 

The variation with true anomaly of the minimum dimensionless derivative 2 of 
thrust-produced velocity impulse with respect to the semimajor axis is pre- 
sented in figure 2(b). The corresponding minimum derivative with respect to 
the period of orbit is obtained by applying the factor 2/3 to the derivatives 
indicated in this figure. It is seen that the minimum local values of these 
derivatives occur at pericenter and the maximum at the point of greatest 
orbital radius. For elliptic orbits the values of either derivative at 
apocenter and pericenter differ by the factor (1 + e)/(l - e) . 

Optimum correction of eccentricity .- The minimum thrust requirements for 
correcting the eccentricity of various orbits are presented in figure 3. From 

2 The dimensionless derivative shown in figure 2(b) has been normalized 
with the circular orbital speed at pericenter, V c . This dimensionless form 

is obtained by multiplying equation (17a) by the factor 

V l f. + 2e cos 0J + e 2 

V7~ = / 1 + e 

C P 

Since V_ is constant for a given orbit and is the same for all orbits 
P 

having the same value of pericenter radius, the orbital position of absolute 
minimum thrust for correction of semimajor axis is obtained from figure 2(b) 
along with a comparison of minimum thrust requirements for various orbits 
having a common pericenter. 
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•x. 


figure 3(a) it is seen that for elliptic orbits the optimum thrust angle for 
increasing eccentricity varies from 0 at pericenter to 180° at apocenter. It 
is noted that the variation of r with Gj shown for e -> 0 is not valid 
for exactly circular orbits for which t e is 0 or 180° at all points on the 
orbit. The optimum thrust direction for parabolic orbits is along the flight 
path; for hyperbolic orbits it varies from along the flight path at pericenter 
to normal to the flight path at infinite radial distance. 

The variation with true anomaly of the minimum dimensionless derivative 
of thrust-produced velocity impulse with respect to the eccentricity is pre- 
sented in figure 3(b). It is seen that the minimum local value of this deriv- 
ative for elliptic and parabolic orbits occurs at pericenter with the maximum 
local value occurring at or near the midpoints of the orbit. For hyperbolic 
orbits the derivative tends toward 0 as the orbital radius tends toward infin- 


ity. In the case of nearly circular orbits the absolute value of 



1.0 at the midpoints of orbit. This derivative has a constant value of 0.5 
at all points for exactly circular orbits. 

Optimum correction of true anomaly . - The minimum thrust requirements for 
correcting the true anomaly (i.e., for slight rotation of the line of apsides 
within the orbital plane) are presented in figure 4. From figure 4(a) it is 
seen that the optimum thrust angle for increase of true anomaly varies from 
-90° at pericenter to 0 at the latus rectum and is rather insensitive to varia- 
tions in eccentricity in this region. In the true anomaly range beyond 90° 
the optimum thrust angle increases to 90° at apocenter for elliptic orbits 
and returns to 0 as the radial distance becomes infinitely large for 
hyperbolic orbits. 

The variation with true anomaly of the minimum dimensionless derivative 
of thrust -produced velocity impulse with respect to true anomaly is presented 
in figure 4(b). For elliptic orbits the minimum value of this derivative 
occurs at or near the midpoint of the orbit with a local maximum occurring at 
both apocenter and pericenter. For parabolic and hyperbolic orbits the 
absolute value of this derivative decreases monotonically with increasing 
absolute value of true anomaly at the thrusting point. 

Optimum correction of orbital radius .- The minimum thrust requirements 
for correcting apocenter or pericenter radius are presented in figure 5. The 
optimum thrust direction for correcting either of these radii (fig. 5(a)) 
varies from the horizontal to the vertical as the point of thrust application 
moves toward the point of orbital radius correction from the opposite point on 
the orbit. Good approximations of the optimum thrust directions in the two 
cases are available along certain portions of the orbit. When rj/r << 1 
the optimum thrust direction for correcting apocenter radius is nearly along 
the flight path, and when rj/r >> 1 the optimum thrust direction for correct- 
ing pericenter radius is nearly^horizontal (cf. appendix D, eqs . (D5) and 
(D6) ) . 

The variation with thrusting-point true anomaly, 0*, of the minimum 
dimensionless derivative of thrust-produced velocity impulse with respect to 
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apocenter or pericenter radius is presented in figure 5(b). In both cases 
this derivative decreases monotonically with increasing angular separation 
between the thrusting and correcting points and attains an absolute minimum 
at the opposite point on the orbit. 

The minimum thrust requirements for correcting orbital radii at the apses 
are special cases of the minimum thrust specifications for correcting any 
orbital radius, r 2 . These general results for correcting orbital radius are 
given by equations (22) and (23) for fixed central angle, 0 2 - 0i, between 
the thrusting and correction points, and by equations (24) and (25) for fixed 
true anomaly, 0 2 , at the correction point. The thrust requirements obtained 
from these expressions are the same when the orbital radius to be corrected is 
at one of the apses. 

It can be shown from equation (22) that the optimum thrust direction for 
correcting orbital radius at fixed central angle from the thrusting point 
varies from the horizontal to a direction normal to the flight path as the 
central angle is reduced from its maximum value to 0. The thrust point loca- 
tion on the orbit for the absolute minimum thrust requirement, however, does 
not generally occur at maximum angular separation from the orbital radius to 
be corrected, as is the case for correcting the apses. This optimum thrust 
point location, in terms of the optimum central angle, 6 2 - 0i*, is presented 
in figure 6 as a function of the true anomaly of the orbital radius to be 
corrected, 0 2 . In general, the point of minimum required thrust along the 
orbit is disposed toward apocenter from the nadir of the correction point 
(i.e., point of maximum angular separation where 0 2 - ©i = 180°). In the 
case of parabolic and hyperbolic orbits the optimum thrust point is located at 
infinite radial distance. 

The abrupt change in location of the minimum thrust point indicated in 
figure 6 at the higher values of eccentricity can be explained in terms of the 
locations of the thrusting and correction points relative to the midpoints of 
orbit. It is found in all cases that the minimum thrust point moves toward 
pericenter as the point of orbital radius correction approaches apocenter. At 
some value of 0 2 the absolute values of 0j* and 0 2 are equal as indicated 
in figure 6. When this equality occurs in the inner part of the orbit as 
delineated by the minor axis (i.e., with ri = r 2 < a) , there is a smooth 
progression of the minimum thrust point toward pericenter as the absolute 
value of 0 2 increases. When this equality occurs in the outer part of the 
orbit (i.e., with r^ = r 2 > a) , the minimum thrust point suddenly shifts 
toward pericenter as the absolute value of 0 2 increases. This effect is 
shown in the curves for the higher values of eccentricity in figure 6. It is 
indicated that simultaneous midpoint locations of the minimum thrust and cor- 
rection points occur for an orbit of eccentricity about 0.7. 

The absolute minimum thrust specifications for correction of orbital 
radius presented in figure 7 were obtained from equations (22) and (23) using 
the optimum values of 0i indicated in figure 6. In figure 7(a) the optimum 
thrust angle results reflect the discontinuous changes in optimum thrust point 
location shown in figure 6. The effects of these changes in the location of 
the optimum thrust point are less apparent in the minimum velocity impulse 
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derivative results shown in figure 7(b). The primary effect is that the 
curves level off rather abruptly at the higher values of eccentricity. 

For correction of orbital radius at fixed true anomaly, the optimum 
thrust angle is 0 when the location of the thrust point is at either apsis. 

From equation (25) it is found that the absolute minimum thrust requirement 
occurs at one or the other apsis, depending on whether the point of orbital 
radius correction is on the inner or the outer half of the orbit. For correct- 
ing orbital radii on the inner half (r 2 < a) , the optimum thrust location is 

at apocenter, and for correcting orbital radii on the outer half (r 2 > a) the 

optimum thrust location is at pericenter. Absolute minimum values of the vel- 
ocity impulse derivative given by equation (25) and based on apocenter or peri- 
center location of the thrust point are indicated by the dashed curves in 
figure 7(b). It is noted that the values of this minimum derivative for cor- 
rection of apocenter or pericenter radius are the same as those presented for 

the case of fixed central angle. 

Correction of elapsed flight time . - Thrust requirements for correcting 
the elapsed flight time to pericenter are presented in figures 8 and 9. The 
results presented in figure 8 are minimum thrust requirements for correcting 
flight time with no constraint on pericenter radius, and those in figure 9 
are thrust requirements for correcting flight time with fixed orbital radius 
at the terminal point. In both cases two different thrust specifications are 
shown for the constraints of fixed central angle, 0 2 - 9i, and fixed true 
anomaly, 0 2 , at the terminal point. 

In figure 8 are shown the variation with thrusting point true anomaly, 

01 , of both the optimum thrust angles, t t and T t ', and the corresponding min- 
imum dimensionless derivatives of thrust -produced velocity impulse with res- 
pect to elapsed flight time to pericenter. It is clear from figure 8(a) that 
the optimum thrust directions in the two cases tend toward and normal to the 
flight-path direction independent of eccentricity as the thrusting point 
approaches pericenter. There is seen to be less variation of the thrust angle 
with eccentricity than of the thrust angle T t '. It is shown that the two 
thrust angles tend toward the same values for the higher values of eccentricity 
as the angular separation and, therefore, the radial distance from pericenter 
increase. From figure 8(b) it is seen that the minimum dimensionless deriva- 
tives tend toward unity and 0 for thrust at angles t t and T t ', respectively, 
as pericenter is approached. As for the optimum thrust angle, there is less 
variation of the minimum derivative with eccentricity for thrust at angle 
than for thrust at angle It is apparent that the two derivatives tend 

toward the same values for the higher values of eccentricity as the distance 
from pericenter increases. 

Somewhat similar effects are to be seen in the results presented in 
figure 9 for correction of the elapsed flight time to pericenter under the 

constraint of fixed pericenter radius. In this case the thrust angles t_ 

r o 

and t r ' are equal and both differ by 90° from the values of x r given in 

figure 5(a). Although the dimensionless derivatives also tend toward unity 
as pericenter is approached under the constraint of fixed central angle, they 
tend toward finite values other than 0 under the constraint of fixed true 
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anomaly at the terminal point. However, as in the results of figure 8 the 
dimensionless derivatives for the two different terminal constraints tend 
toward the same values as the angular separation from pericenter increases. 


Application of Minimum-Thrust Results to Interplanetary Guidance 

The guidance problems considered here occur during interplanetary 
missions when the spacecraft has not been precisely injected into the planned 
transfer orbit, and subsequent flight-path corrections are needed. It is con- 
venient to separate guidance for such missions into a midcourse phase and an 
approach phase depending on whether the vehicle is outside or inside the 
sphere of influence of the target planet. When the guidance is based on the 
use of a reference trajectory, it is termed implicit guidance in contrast to 
explicit guidance for which no reference trajectory is required. 

During the midcourse phase flight-path corrections are made to ensure 
that the spacecraft is on a collision course with the target planet or a point 
slightly removed from the planet. Both fixed and variable time of arrival 
guidance have been considered for this purpose (e.g., see refs. 8-14). In the 
case of implicit guidance there is little difference in the guidance computa- 
tions whether the time of planet intercept is considered fixed or variable. 

In either case it is necessary to know the flight-path deviations from the 
reference trajectory at or near the terminal point. 

During the approach phase final corrections of the vehicle flight path 
are made prior to actual interception of the target planet. Along with cor- 
rection of the orbital plane, control of the time and/or location of peri- 
center in a planetocentric frame of reference will usually be required. For 
this purpose a form of fixed radius of pericenter guidance, based on the in- 
plane pericenter deviations from a reference trajectory, may be used. As is 
the case whenever impulsive thrust corrections are applied, three elements or 
parameters of the orbit can be corrected with each application of thrust. 

The manner in which the minimum thrust results of the present analysis 
apply to several specific types of implicit guidance during the midcourse and 
approach phases of interplanetary flight is next described. In both phases 
only the in-plane components of the total required thrust correction are con- 
sidered, since the out-of-plane corrections controlling the plane of orbit are 
generally independent of the in-plane corrections and can be treated sepa- 
rately. It is assumed that navigation information is available as to the in- 
plane flight-path deviations of the uncorrected spacecraft orbit from the 
reference trajectory at the terminal point. A summary of the guidance 
computations outlined in this section is presented in table II. 

Midcourse guidance .- During the midcourse phase of interplanetary flight 
the spacecraft is between the spheres of influence of the launch and target 
planets. For this phase of flight the spacecraft motion relative to the Sun 
is essentially two-body planar motion so that the first-order minimum thrust 
relationships of the present analysis are suitable in connection with linear 
guidance about a reference trajectory. This guidance may be such as to result 
in interception of the target planet at the reference intercept time (fixed 
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time of arrival) or at some time before or after the reference intercept time 
(variable time of arrival) depending on the particular constraints to be 
satisfied. Both cases are considered below. 


The general situation prior to midcourse guidance action is illustrated 
in figure 10 which shows the orbits of the spacecraft and the planet in the 
vicinity of the planned intercept point. Interception of the target planet 
is designed to occur at time t 2 with the planet at position P 2 . The uncor- 
rected flight path of the spacecraft is radially displaced a distance fir 2 
from the design intercept point and crosses the radial line through this point 
at a time different by fit 2 from the design value. In fixed time of arrival 
guidance the in-plane flight path corrections required are precisely these 
deviations from the reference trajectory at the terminal point. The direction 
and magnitude of the in-plane velocity impulse required to correct the flight 
path, in terms of the end-point deviations, are obtained from equations (34) 
and (35) expressed as 


A V - AV 


cos (t - x ) AV, 

T v t T J 1 


tan(x - x r ) - av t sin(x t - x r ) 


AV 


+ tan(x t - x ) (36) 

L L o 


AV = AV t /cos(x - x„) 

T L T 1 


(37) 


where, since d(t 2 - t|) 
dV T 


AV, 


dr 2 


fir 


2 > 


= dt- 


AV, 


dV 


T t 


dt; 


fit 


2 > 


AV, 


dV, 


dtc 


fit; 


The optimum thrust angles T r and x t are given in the text by equations (22) 
and (26) and the minimum derivatives dV T ^/dr 2 and dV Tt /dt 2 by equations (23) 

and (27). The derivative dV T /dt 2 is given by equation (29a) with the 

r o 

thrust angle x„ different from x r by 90°. In all cases the thrust angles 
L o 1 

and derivatives are evaluated on the reference trajectory. 


In variable time of arrival guidance about a reference orbit, intercep- 
tion of the target planet can occur at any point along its orbit in the vicin- 
ity of the nominal intercept point. Such a point is represented in figure 10 
by P 2 ' separated from the nominal intercept point, P 2 , by the central angle 
fi9 2 . The in-plane components of orbit deviation from planet interception at 
point P 2 ' can be expressed in terms of the components of deviation from 
interception at point P 2 approximately as 

fir 2 ' = fir 2 - r 2 (tan y s - tan 1 


fit 2 ' = fit 2 


u b 



( 38 ) 
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where y s , u s and y b , i% are nominal values of flight-path angle and hori- 
zontal component of velocity for the spacecraft and planet, respectively, at 
the nominal intercept time, t 2 • With substitution of these adjusted values of 
6r 2 and 6t 2 in the thrust vector computation above, the required in-plane 
corrective velocity impulse is obtained for interception of the target planet 
at point P 2 '. It is assumed in this procedure that the optimum thrust angles 
and minimum derivatives in the computation are invariant, to first order, for 
small values of S0 2 . The computation reduces to the fixed-time-of-arrival 
result when 66 2 is 0. 

For given values of the components of deviation, 6r 2 and 6t 2 , from planet 
interception at the nominal or design intercept point, P 2 , there is a partic- 
ular value of 66 2 at each point along the uncorrected vehicle flight path 
for which the magnitude of the in-plane corrective velocity impulse for planet 
interception is a minimum. This value of 60 2 is obtained, after introduc- 
tion of equation (38) into equation (37), by setting the derivative of AV X 
with respect to 60 2 equal to 0 and solving for 60 2 . It is found that 

dV T dV x 

[D r - cos(x r - T.j.) ] 6r 2 + [D-j- - D^. cos(Xj. - x-j-)] — 6t 2 

60 2 = , _ , — — 

D r 2 + D t 2 - 2D r D t cos(x r - x t ) 

(39) 


where 

dV T 

D r = r 2 (tan y g - tan y b ) ^ 

Dt . (i . V) p, 

1 U b V U S / dt 2 

The in-plane thrust requirements (magnitude and direction of velocity impulse) 
for the minimum- thrust maneuver are obtained by simply inserting this value 
of 60 2 into the expressions for 6r 2 * and 6t 2 f and substituting these 
adjusted terminal point deviations for 6r 2 and 6t 2 in equations (36) and 
(37). 

Approach guidance .- During the approach phase of interplanetary flight 
the spacecraft is within the sphere of influence of the target planet. Here 
the vehicle motion relative to the planet is adequately defined by the two- 
body equations of motion, and minimum thrust relationships derived from these 
two-body equations can be used in connection with linear guidance about a 
reference trajectory. The type of implicit guidance obtained is essentially 
the same as that given for the midcourse phase except that additional 
simplification is possible. 

Various forms of implicit guidance have been considered in the literature 
for the approach phase. These include fixed time of arrival guidance wherein 


29 



the spacecraft is constrained to pass through the reference terminal point 
(usually pericenter of the reference trajectory) at the design arrival time. 

In this case corrective thrust computations are identical to those already 
described for this form of guidance in the midcourse phase. Although the 
pericenter radius is not explicitly controlled with this type of guidance, the 
deviation of pericenter radius from its reference value will usually be small, 
and the time of pericenter passage will not differ greatly from the reference 
arrival time. 

When explicit control of the pericenter radius is desired, some form of 
fixed radius of pericenter guidance is used. Along with the constraint on 
pericenter radius a second constraint must be imposed to uniquely define the 
in-plane corrective thrust vector. Three such secondary constraints which 
merit consideration are: (1) minimum corrective thrust, (2) fixed time of 

pericenter passage, and (3) fixed location of pericenter. The computation of 
in-plane corrective velocity impulse vectors for approach guidance with these 
constraints is indicated below. 


The general situation prior to approach guidance action is illustrated in 
figure 11 which shows the in-plane deviation of the uncorrected flight path of 
the spacecraft from the reference trajectory. Pericenter of the uncori’ected 
orbit has components of deviation fir D and 60 p from the reference location 
along with an error, St-,, in the time^of pericenter passage. From equation (8) 
corrections to these pericenter parameters due to a velocity impulse at 
angle t of magnitude AV X can be expressed as 




6r p = 


66 p - 


fit = 
P 


AV t cos(t - T r ') 

d V'' dr P P 

r p 

AV t cos (t - Tg) 

dv X0 /de 1 

AV X cos (t - t^) AV X cos (t 




(40) 


dV 


T t; /dt e 


dV 
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T » 

r o 


dt 


COS(T r I - T t t) 


J 


where the optimum thrust angles, x r ' , T t ', and t , are given by equations (D6) , 

(E3) , and (20), respectively, and t^e minimum derivatives dV x i/dr , 

r P P 

dV x ,/dt , and dV x /dS} by equations (DIO), (E5) , and (21). In the alternate 


expression for fit the thrust angle T r 


is different from r r ' by 90°, 


L o ■ L p 

and the derivative dV T ,/dt is obtained from equation (E8) . 

r o P 

It is possible to correct any pair of the three pericenter parameters 
expressed in equation (40) with a single application of impulsive thrust. The 
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thrust requirements for such a correction are obtained by solving simultan- 
eously the pertinent pair of equations for t and AV t with the result given 
in the form of equations (34) and (35). In the following three cases of 
approach guidance considered, one of the two parameters being corrected in 
each case is the pericenter radius. The other parameter being corrected or 
the secondary constraint imposed in each instance is one of the three listed 
above. 

For approach guidance which corrects the pericenter radius deviation with 
minimum thrust, the direction and magnitude of the in-plane corrective veloc- 
ity impulse are given by 


T 


AV 


T 




6r P 


In this case the time of pericenter passage and the angular position of peri- 
center after thrusting will generally differ from the reference values. These 
differences are readily computed by introducing the above values of t and 
AV t into equations (40) and adding the resulting changes in tp and 0p to 
the initial errors. 

For approach guidance which corrects the deviations in pericenter radius 
and time of pericenter passage, the direction and magnitude of the in-plane 
corrective velocity impulse are given by 


AV 


V " AV ^r' COs(T tp - T r'> dV ^r! dr p 6t r 


tan(x - T r t) = 


AV X i sin (x t i - x r 0 


dt T 


dV x , Sr 

T r p P 


■cot(x t ^ 


T r,0 


AV 


AV x r ’ 

P 

cos fx - X i) 
r ' 


dV T ' 


6r 


dr cos(t - t ») 


In this case the angular position of pericenter after thrusting will generally 
differ from the reference value. This difference is readily computed using 
equations (40) and the initial error in orbit orientation. 

For approach guidance which corrects the location of pericenter in the 
plane of orbit, the direction and magnitude of the in-plane corrective 
velocity impulse are given by 


31 


A y - AV T , cos (Tg - T r ,) dV dr J9 


tanfx - t ' ) 


AV x . sin(T 0 - t •) d 0 j dV T « 6 r 

x -n P j-n r 


— xzr- -cot(T c 


T rP 


AV 


T r 1 


dV , 

L -r 


6r 


P P ”p 

AV = 7 TV = -j 7 r 

T COS (T - T 'J dr cos (x - T 0 

P T 


In this case the time of pericenter passage will generally differ from the 
reference value. This difference is readily computed using equations (40) and 
the initial error in t . 


SUMMARY OF RESULTS 


From a first-order analysis of the Keplerian-Newtonian equations of 
motion, expressions have been developed that give the minimum thrust require- 
ments for effecting small changes in certain elements or parameters of two- 
body orbits. The thrust requirements are presented in terms of the optimum 
thrust angle and the derivative of the thrust-produced velocity impulse with 
respect to the orbital parameter to be corrected. A summary of such thrust 
specifications for single parameter correction is presented in table I. 

The optimum thrust angles for correction of various orbital parameters 
are invariant, to first order, with the amount of thrust applied. For correct- 
ing either the semimajor axis or the period of orbit, the optimum thrust direc- 
tion is along the flight path at all points along the orbit for all values of 
eccentricity. For correcting the eccentricity of elliptic orbits, the optimum 
thrust direction is along the flight path at apocenter and pericenter. At 
other points along elliptic orbits optimum-angle thrust for correction of 
eccentricity tends toward the flight-path direction as the value of eccentric- 
ity increases to 1 (parabolic orbit). A minimum amount of optimum-angle 
thrust is required at pericenter for correcting either the eccentricity or the 
orbital period. In the case of correcting true anomaly or rotating the line 
of apsides of an orbit, the absolute value of the optimum thrust angle varies 
from 90° at apocenter and pericenter to 0 on the semilatus rectum. For small 
amounts of rotation, a minimum amount of optimum-angle thrust is required near 
the midpoints of the orbit. To first order, the optimum thrust angle for mini- 
mizing the change in any parameter is different by 90° from the optimum thrust 
angle for maximizing the change in that parameter. 

The optimum thrust direction for maximum change in the orbital radius at 
a given point along an orbit varies from the horizontal at the nadir point on 
the orbit to a direction normal to the flight path at the given point. The 
amount of optimum-angle thrust required in correcting the orbital radius at a 
given point along the orbit increases from a minimum at some point on the 
opposite side of the orbit (which is generally disposed toward apocenter from 
the nadir point) to a maximum (infinite value) at the point of desired radius 
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change. The thrust direction required in correcting the elapsed flight time 
to a given point along an orbit is normal to the thrust direction required in 
maximizing the orbital radius change at this point. The thrust magnitude 
required in correcting the elapsed flight time to a given point (as a percent 
of the value of orbital velocity at the given point) tends toward the 
percentage change in elapsed flight time as the given point is approached. 

In general, it is possible to correct two independent parameters of a 
given orbit with a single application of impulsive thrust. The thrust 
requirements for such a correction can be formulated in terms of the minimum 
thrust requirements for correcting each parameter separately. These results 
are found to have application to the problem of spacecraft guidance during 
interplanetary flight. The application has specific relevance to implicit 
guidance within the plane of orbit and is essentially the same whether the 
spacecraft is in the intermediate stage of its flight or is in the final stage 
within the sphere of influence of the target planet. In either case naviga- 
tion information is required of the in-plane deviations of the projected 
vehicle flight path from a design or reference trajectory at the terminal 
point. For midcourse guidance these terminal-point deviations consist of the 
differences in orbital radius and elapsed flight time from reference values 
at the desired point of planet interception. This desired point of intercep- 
tion can be any point along the orbit of the planet in the vicinity of the 
nominal or design intercept point. For approach guidance, after penetration 
of the sphere of influence of the target planet, the terminal-point deviations 
consist of the differences from reference values in pericenter location or in 
pericenter radius and time of pericenter passage. A summary of in-plane 
guidance computations based on use of the minimum- thrust results of the 
present analysis is given in table II. 
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APPENDIX A 


SUMMARY OF USEFUL KEPLERIAN-NEWTONIAN RELATIONS FOR TWO- BODY MOTION 
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APPENDIX B 


ALTERNATE FORMULAS FOR CORRECTING ORBITAL RADIUS 


The following expressions are alternates to equations ( 13 ) and ( 22 ) 
through ( 25 ) for thrust requirements for correcting orbital radius r 2 . These 
results are obtained using Vi and as dependent variables defining the 
spacecraft velocity vector at the thrusting point with the thrust angle, 
t - Yi, specified relative to the flight path. The results consist of thrust 
specifications for both arbitrary and optimum angles of thrust. The general 
expression of these results is given by equations ( 4 a), ( 5 a), and (6a) when 
the terminal point constraint on 0 2 is specified (either d0 2 = d0j or 
d0 2 = 0). 

For correcting orbital radius at a fixed central angle, 0 2 - 0 X , 
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For correcting orbital radius at fixed true anomaly, 0 2 , 


v l dr 2 

^2 dV^ 


2(1 - cos 0j cos 0 2 ) 
1 + e cos 0 2 


cos (t ' 


Yi) 


[ (1 + e 2 )cos 0 2 + 2e]sin 0]^ 
(1 + e cos 0 2 ) (1 + e cos 0j) 


sin(x ' 


Yi) 


(B 4 ) 


37 



[(1 + e 2 )cos 0 2 + 2e]sin ©i 
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APPENDIX C 


ALTERNATE FORMULAS FOR CORRECTING ELAPSED FLIGHT TIME 


The following expressions are alternates to equations (15) and (26) 
through (29) for correcting elapsed flight time t 2 - tj. As in appendix B 
the results are obtained using V* and yi as dependent variables defining the 
spacecraft velocity vector at the thrusting point with the thrust angle, 
t - yi» specified relative to the flight path. Again the results consist of 
thrust specifications for both arbitrary and optimum angles of thrust as given 
in general form by equations (4a), (5a), and (6a) when the terminal point 
constraint on 02 is specified (either d 02 = d0i or d ©2 = 0). 


Vi d(t2 - tj) 
t2 - ti d\T 


' / P 3 /v> 

t2 . tl t Y l cos (t - yi) + Y 2 sin(x - Yj)] (Cl) 


where 


Y 1 

y 2 


2B 

e 


J - - — sin 0! 


A sin 0i - B 


1 + e 2 

2 + cos 0i 


1 + e cos 0i 


A 


(2 + e cos 02 )sin ©2 (2 + e cos 0i)sin 0i 

(1 + e cos © 2 ) 2 (1 + e cos ©i) 2 


B = (1 + e .cos 02 ) 2 - (1 + e cos 0i) 2 


3(1 + 2e cos 0i + e 2 ) t 2 - ti 2A(e + cos 0i) 

j = — 

1 - e 2 /p 3 / p 1 - e 2 


In the case of parabolic orbits J is replaced by its limiting value defined 
by 


lim J = 2L(1 + cos 0i) 
e->-l 

where L is given in the text with equation (15) . 
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The minimum thrust requirements for correcting elapsed flight time to 
traverse a fixed central angle with optimum-angle thrust are given by 


Y 2 Ae sin 0i - B [2e + (1 + e 2 )cos 0 ^ ] 
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For correcting elapsed flight time to arrive at a fixed true anomaly, 
0 2 , B is replaced in the above expressions by B' where 


B' = -(1 + e cos 0i)“ 2 


As is indicated in the text, thrust at angle or t~ ' will correct 
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the elapsed flight time to a given terminal point on the orbit without chang- 
ing the orbital radius at the terminal point to first order. In the one case 
(x r ^) the terminal point is at a fixed central angle, 0 2 - 0^, while in the 

other case (t ') the terminal point is at a fixed true anomaly, 0 2 . When 
r o 

these thrust angles are inserted in equation (Cl) , it is found that 
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where A, B, B', and J are defined above. 
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APPENDIX D 


FORMULAS FOR CORRECTING APOCENTER OR PERICENTER RADIUS 


The following expressions are special cases (02 = tt, 0) of the general 
thrust requirements for orbital radius correction presented in the text and in 
appendix B. These thrust specifications for correction of apocenter and peri- 
center radii are independent of the terminal point constraint on 0 2 due to 
the insensitivity of orbital radius at 02 = it or 0 to small changes in 02 - 


For thrust at arbitrary angle there are obtained from equation (13) 
V dr 2(1 + cos 0i) - e sin 2 0i 


r dV 
a t 


1 + e cos 0 i 


cos x + sin 0 i sin t 


(Dl) 


V dr_ 2(1 - cos 9i) + e sin 2 0j 


r dV 
P T 

and from equation (Bl) 


1 + e cos 0i 


cos t - sin 0 i sin t 


V i dr 2(1 + cos 0 i ) 
a 


r a dV T 


cos (t - yi) + 


(1 - e)sin 0i 
1 + e cos 0i 


sin (t - Yj) 


(D2) 


(D3) 


V i drp 2(1 - cos 9 i ) 
r~ dV = 

p T 


1 + e 


(1 + e)sin 0i 

cos(T " Yl) - i% e 'cos TT 5in(T “ Yl 5 


(D4) 


Optimum thrust angles for correcting apocenter or pericenter radius are 
obtained from equation (24) as 
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Minimum dimensionless derivatives of thrust -produced velocity impulse 
with respect to apocenter and pericenter radii are obtained from equation (25b) 
as 
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APPENDIX E 


FORMULAS FOR CORRECTING FLIGHT TIME TO APOCENTER OR PERICENTER 


The following expressions are special cases (0 2 = tt, 0) of the general 
thrust requirements for elapsed flight time correction presented in the text 
and in appendix C. In obtaining these expressions the constraint of fixed 
true anomaly at the terminal point, 0 2 , has been imposed. For the sake of 
brevity, in each instance a single formula is presented that is valid for both 
cases of terminal point location. It is understood that the subscript 2 is 
to be replaced by the subscript a when the terminal point is apocenter and 
by the subscript p when the terminal point is peri center. 


For thrust at arbitrary angle there is obtained from equation (15) 
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and from equation (Cl) 
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Optimum thrust angles for correcting elapsed flight time to apocenter or 
pericenter are obtained from equations (26) and (C2) as 
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Minimum dimensionless derivatives of thrust-produced velocity impulse 
with respect to elapsed flight time to apocenter or pericenter are obtained 
from equations (27) and (C3) as 
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t 2 ~ T t^ 


"I 2 pfl+sin 2 eO-cos 9i / 3 /„ „ 1 : 

l) + : — T + 3 e 2 sin ®l 

1 J L 1+e cos 6i t 2 -ti ‘J 

e(t 2 - ti)/>^) 3 /ij 


(E5) 


Vx d(t 2 - 1 1 ) 


1 + 2 e cos 6 x + e 2 "! 2 !" 2 sin 6x 3 e(t 2 - t x 3 ”| 2 cos 2 0 X 


1 - e 2 


1+e cos 0 x 


/p 3 /y 


Cl + e cos 0 1 ) 2 

(E6) 


The corresponding dimensionless derivatives for thrust in the direction 
that, to first order, does not change the orbital radius at the terminal point 
(t = T r *) are obtained from equation (29b) as 


t - ti 
a 1 


dV T , 


V a d(t a ' t l ) 


+ 1-e /p 3 /u 


cos 0 i 3 - e sin 2 0 i ] 2 + (1 + e cos 0x) 2 sin 2 0x 


2 /l + cos 6 1 - e sin 2 6x \ 

e \ 1 e cos 0^ / 


*a - tl 

+ 3(1 - e) — 7 sin 6x 

/p 3 /y 


(E7) 




dV, 


1 - e "p 


~ ti 


Tl o = + 1 + e /p 3 /y 


V n d(t p " t l ) 


/[2(1 - cos 0 x ) + e sin 2 0 1 J 2 + (1 + e cos 0]) 2 sin 2 0j 


1 - cos 0 x + e sin 2 0x \ t p - tx 


1+e cos ©x 


0 


+ 3(1 + e) 


/p 3 /u 


sin 9 j 


CE8) 
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f 


and from equation (C5) as 


*a - ti 

Vi 


dV T , 
r o 




/ A (1 + cos 0i) 2 (1 + e cos 6 1 ) 2 + (1 - e) 4 sin 2 0! 


d(t a - to 


1 + 2e cos 


cos e l + e [2 /I + cos 
1 + e Je \ 1 + 


1 + cos 0] - e sin 2 0i 


e cos 0i 


\ t a - t! ] 

-j + 3(1 - e) == sin 0j 
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(E9) 


*P - tl 
Vi 


dV T 
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= ± ... 

1 + 2e cos 0i + e 2 p 2 /I - cos ©i + 

1 - e |_e V 1 + ec 
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-V 
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TABLE I.- THRUST SPECIFICATIONS FOR CORRECTION OF ORBITAL 

PARAMETERS SINGLY 


Semimajor axis, a, 
or 

period of orbit, T 


Constraint 


Minimum 

thrust 


Eccentricity, e Minimum 

thrust 


True anomaly, Qj Minimum 

thrust 


Thrust angle Velocity impulse derivative 

Function I Equation Function Equation 


Minimum thrust 
point on orbit 


(16a) 

(16b) 

a dV T a T dV r a 
Vj da ’ V! dT 

a dV r a T Ha 
V p da ’ V p dT 

(17a) 

(17b) 

Pericenter 

(18a) 

(18b) 

1 dV -e 
V : de 

i dv x 
1 T e 

Uj de 

(19a) 

(19b) 

Peri center 
(for e < 1) 

Infinity 
(for e > 1) 

(20a) 

i dv ^e 
V, de) 

(21a) 

Near midpoint 
(for e < 1) 

(20b) 

Ui d0 i 

(21b) 

Infinity 
(for e > 1) 


I 

I 


Orbital radius, r 2 Minimum 

thrust , 

9 2 - 9i = const. 


tan(x r -Yi) (B2) 

tan x (22) 


180° > 0! > 0 2 + 180 
(for 0 2 0) 

-180° < 0j < 0 2 - 18 

(for e 2 > 0) 


Minimum 

thrust, 

0 2 = const. 


tan(x »-yi) (B5) 


*•2 dv T ; 

V! dr 2 


Apocenter 
(for r 2 a) 

Peri center 
(for r 2 21 a ) 


Elapsed flight 
time, t 2 - ti 


Minimum 

thrust, 

0 2 - Bi - const. 


Minimum 

thrust, 

0 2 = const. 


r 2 = const., 

0 2 - 0! = const. 


r 2 = const. , 
0 2 = const. 





t 2 - 

tl 

dv x t 

tan(T t - 

Yi) 

(C2) 

Vx 

— 

d(t 2 - ti) 




t 2 - 

tl 


tan x^. 


(26) 

U 1 


d(t 2 - ti) 




t 2 - 

tl 

dv T . 

tan(x t - 

1 1 'Y l) 

1 (C2) i 

B = B ' 

1 


d(t 2 - t]) 

tan x t 


(26) , 
B = B' 

t 2 - 
ui 

tl 

dV T^ 

d(t 2 - ti) 


tan(Tr 0 -Yj) 
tan t 

r o 

(28a) 

t 2 - t, dV Tr 0 

Vj d(t 2 - ti) 

t 2 -t, av Tro 

Ui d(t 2 - ti) 

(C4) 

(29a) 

tan(T r ^- Yl ) 

l 

t 2 - t! dV X r ^ 

V) • d(t 2 - ti) 

(C5) 

tan x_ ’ 
0 

(28b) 

ta-ti dv ^r 0 

Ui d(t 2 - tl) 

(29b) 
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TABLE II.- APPLICATION OF MINIMUM THRUST RESULTS TO IN-PLANE INTERPLANETARY GUIDANCE 


Type 

of 

guidance 


Guidance 

objectives 

or 

trajectory 

constraints 


Orbital parameter 
corrections required 


Direction and magnitude 
of 

corrective velocity impulse 


Source of quantities used in guidance 
computations 


Optimum thrust angle and 
minimum velocity impulse for| 
single parameter correction 


Minimum 

thrust 

produced 

velocity 

impulse 

derivative 


Fi xed 
time of 
arrival 


Planet 

interception 
at time t 2 


Terminal point 
deviations, 6 T 2 
and 6 t 2 


r av /av 

T “ T + tan 1 — : — ^r— 

l_sin(-r t - T r ) 


- COt (t t - T r )J 
aV T 

av = c- 


t COs(t - T ) 


T r from 
eq. ( 22 ) 


T t from 
eq. (26) 




■ from 


dr 2 
eq. (23) 


dV T 


- 6t 2 


dt 2 
eq. (27) 


- from 


Variable 
time of 
arrival 


Planet 

interception 
at time 
near t 2 


6 r 2 ' = 6 r 2 


-r 2 (tan y g - tan Y b ) 68 2 


6 t 2 ' 


6t 2 

-II (l- 02) 
u b V V 


dV T 

AV = -r~ 6 t 2 
T r dr 2 


4 \ = dtr iV 


Same 


60 2 


Planet 

interception 
with minimum 
thrust 


Same as above , 
60 2 from eq. (39) 


above 


above 


above 


Fixed 

radius 

of 

I pericenter 


Peri center 
radius 

control with 
minimum 
thrust 


Fixed values 
of peri center 
radius and 
time of 
pericenter 
passage 


Fixed 

location of 
pericenter 


Pericenter radius 
deviation. 6 r 


T r ' from 
P 

eq. (D 6 ) 


dV - 

~~dr~ 


dV T 


eq. (DIO) 
or eq. (D12) 


Pericenter 
deviations , 
6 r p and 6 t p 


W T ./AV 

lP.- r P. - 

LsinCiy - t r 0 


- cotCx’ - ) 

P PJ 


Peri center 
deviations, 
5 r and 60„ 


P PJ 

COS (t - T ') 

p 


x r 1 from 
P 

eq. (D 6 ) 


* ' from 
r P 

eq. (E3) 


r * V T _/AV t - 
9 T P 

L sin(i 0 - t 0 

n p 

- COt (T e - T r ’)j 


t = r ' + tan 
r P 


AV r 


T COS(T - T 1) 

r P 


. 1 from 
eq P (D 6 ) 


t a from 
eq. ( 20 ) 


% % 


dv s 

TT* 4 *, 

p F 


dV T . 


- from 


av t . 


r p 


dV T . 
r P 

-IF- 2 ' 


dv T 

AV = *■?■? — — 68 

T e de 1 p 


dr p 
eq. (DIO) 
or eq. (D12) 

dV T , 

-rr^ from 

% 

eq. (ES) 
or eq. (E 6 ) 


d S 

-jj—from 

eq. (DIO) 
or eq. (D12) 

dV T 


01 


■ from 


eq. (21) 


J 





(a) Parameters of elliptic orbits. 

Figure 1.- Convention of notation used in analysis. 
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(b) Minimum derivative of thrust -produced velocity 


Figure 3.- 



impulse with respect to eccentricity 
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(b) Minimum derivative of thrust -produced velocity impul 

Figure 4.- Concluded. 



with respect to true anomaly. 
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Thrusting point true anomaly, 5 |f deg 

(a) Optimum thrust angles, T r and T r . 

a p 

Figure 5.- Minimum thrust requirements for correction of apocenter and pericenter radii. 
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(b) Minimum derivative of thrust-produced velocity impulse wit 

radius . 


Figure 5.- Concluded. 



respect to apocenter or pericenter 




Figure 6.- Optimum location in orbit for correction of orbital radius under 
constraint of fixed central angle, 62 - 6 1 . 
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Correction point true anomaly, 0 2 , deg 

(b) Absolute minimum derivative of thrust-produced velocity impulse with respect to orbital radius at 

any point along orbit. 


Figure 7.- Concluded. 





Thrusting point true anomaly, Q it deg 


(a) Optimum thrust angles, x and x t ». 

Figure 8.- Minimum thrust requirements for correction of elapsed flight time 

to pericenter. 
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Thrusting point true anomaly, a () deg 

(b) Minimum derivative of thrust-produced velocity impulse with respect to elapsed flight time to 

pericenter. 


Figure 8.- Concluded. 
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Figure 9.- Dimensionless derivative of thrust -produced velocity impulse with respect to elapsed flight 
time to pericenter under constraint of fixed value of pericenter radius. 
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Figure 11.- General situation prior to approach guidance action showing 
in-plane deviation of uncorrected spacecraft flight path from reference 
trajectory. 
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